Effective medium theory of semiflexible filamentous networks 
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We develop an effective medium approach to the mechanics of disordered, semiflexible polymer 
networks and study the response of such networks to uniform and nonuniform strain. We identify 
distinct elastic regimes in which the contributions of either filament bending or stretching to the 
macroscopic modulus vanish. We also show that our effective medium theory predicts a crossover 
between affine and non-affine strain, consistent with both prior numerical studies and scaling theory. 
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Semiflexible polymer networks form a distinct class of 
gels whose mechanical properties are important for both 
biophysical and materials research. These cross-linked 
polymer networks differ substantially from the flexible 
polymer gels and rubbers [l| due to the rigidity of the in- 
dividual polymers ^ . Because the thermal persistence 
length of the constituent filaments is much longer than 
the typical distance between cross-links, these materi- 
als can store elastic strain energy in both stretching and 
bending deformations of the filaments. The cytoskele- 
ton, a complex assembly that includes stiff filamentous 
proteins present in most eukaryotic cells, is an especially 
common example of such a semiflexible network [4| . Such 
networks dominate the mechanical properties of the cy- 
tosol and are at the heart of cellular force production and 
morphological control. 

Theoretical studies of the elastic response of randomly 
cross-linked, stiff, filamentous networks have recently un- 
covered a surprising cross-over between distinct mechan- 
ical regimes of these semiflexible networks Q. For 
given filament elastic parameters there is a transition 
from strain energy storage in filament stretching modes 
at higher network densities to filament bending modes 
in more sparse networks. This transition is accompanied 
by a change in the geometry of the deformation field over 
mesoscopic lengths. At higher densities the network de- 
forms affinely as expected from continuum elasticity the- 
ory while at lower densities, where the elastic energy is 
stored in bending modes, the network deformation field 
is nonaffine over mesoscopic distances. Recent experi- 
ments 0, @ support the existence of this affine (A) to 
nonaffine (NA) cross-over. However, a fundamental un- 
derstanding of the relation of the network architecture 
and individual filament mechanics to the collective elas- 
ticity of the network remains elusive and prior theoretical 
work has been primarily numerical. 

In this letter we develop an analytical model of the me- 
chanical response of two-dimensional, disordered, semi- 
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flexible networks. We introduce a mechanical mean-field 
or effective medium theory of the system that allows us 
to calculate the elastic response of the system to uni- 
formly imposed as well as wavenumber-dependent strain 
fields. From this mechanical response we identify an 
A/NA cross-over and obtain a phase diagram of the sys- 
tem showing the regimes of affine and non-affine behav- 
ior. Our study also demonstrates the presence of a nat- 
ural length scale controlling the A/NA cross-over that 
corresponds well with prior results from simulation and 
scaling theory 0]. 

We study a model two-dimensional system constructed 
as follows. We arrange infinitely long filaments in the 
plane of a two-dimensional hexagonal lattice so that at 
each lattice point three filaments cross. In this way each 
lattice point is connected to its nearest neighbor by a sin- 
gle filament. A sketch of the network is shown in Fig. [1] 
The filaments are given an extensional spring constant a 
and a bending modulus k. The cross-links at each lat- 
tice site do not constrain the angle between the crossing 
filaments. We introduce finite filament length L into the 
system by cutting bonds with probability 1 — p, where 
< p < 1, with no spatial correlations between these 
cutting points. This generates a broad distribution of 
filament lengths with a finite average (L), while intro- 
ducing quenched disorder. We then study the mechanical 
response of this disordered network in the linear response 
regime. The main assumption in our theory is that the 
depleted network has the same mechanical response as a 
uniform network with effective elastic constants am and 
Km. These are determined by requiring that strain fiuc- 
tuations produced in the original, ordered network by 
randomly cutting the filaments have zero average. Here, 
we do not explicitly consider thermal fluctuations, whose 
role in determining the longitudinal compliance of fila- 
ments has been discussed before 9]. These thermal ef- 
fects can be incorporated in the present model through 
a renormalized parameter a 

The elastic energy of the strained network, arising from 
bending and stretching of the constituent filaments, can 
be written in terms of the displacement vector u, at each 
lattice site i. To quadratic order in u the stretching {Eg) 
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FIG. 1; (color online) Schematic figure of the filament net- 
work. The solid red lines represent the undeformed fila- 
ment network, while the dashed blue lines show the deforma- 
tion field having wavevector q and displacement amplitude 
u (shown in the upper left corner of the figure). The black 
arrows show the displacement field at each lattice point. This 
perfect lattice is disordered by making randomly placed cuts 
in the infinitely long filaments. These are not shown. 



and bending (Ef,) energies are 

Eb = :;nR^'^ ^ {uih X - Uij X rihY 



(1) 
(2) 



where R is the lattice constant, is a unit vector di- 
rected from the to the j*^ equilibrium lattice site, 
and Uij is the difference in the strain field between those 
lattice sites. 

It is now simple to determine the collective elastic 
properties of the perfect lattice; doing so for the dis- 
ordered lattice generated by randomly cutting the fila- 
ments is less trivial. We determine the spring constant 
and bending modulus of a spatially uniform effective sys- 
tem [l^ that reproduces the mechanics of our disordered 
system in an average sense as described below. 

We first apply a uniform dilation to the uniform system 
with spring constant am so that all bonds are stretched 
by d£. There is no bending deformation. If we now re- 
place a single filament segment connecting points (say) 
i and j (see Fig. [1]) by one of spring constant a', the 
virtual force needed to the fix positions of i and j is 
/ = d£{am — oi'). If / were applied to the same segment 
in the unstrained network the resulting change in length 
would be J /{am/a* ~ a„i + a'), where a* (0 < a* < 1) 
is a network material parameter that includes the con- 
tribution of the elasticity of the entire network. It may 
be written in term of the dynamical matrix of the lattice 
D{q) as 



«* = ^E^^ [Ds{q)-D-\q)\ 



(3) 



where the sum is over the first Brillouin zone. Here 
D[q) = Ds{q) + Diy{q), where Dg^tiQ) define the stretch- 
ing and bending contributions, respectively, to the full 
dynamical matrix and are given by: 



(4) 



Db{q) = K,nR ^^[4(l-cos(q.fy)) 

-(1 - cos(2q.ry ))] (J - rij-r^- ) (5) 

with I the unit tensor and the sums are over nearest 
neighbors Note that for small q, Dt ^ q^ and 

Ds ~ q^ have the expected wavenumber dependencies 
for bending and stretching. 

From linearity it follows that the extra displacement 
5u of the segment ij due to the change in that filament 
segment's spring constant in the dilated network is the 
same as its extension in response to the force / being 
applied to it. Therefore this additional displacement or 
'fluctuation' is: 



5u ~ 



OLra/a* 



a' 



(6) 



We now average this extra displacement over the ensem- 
ble of possible filament substitutions, with the statistical 
distribution of longitudinal spring constants: 



P(a') = pj(a - a') + (I - p)(5(a'), 



(7) 



where 1 — p is the probability of a cut bond and (5(. . .) is 
the Dirac delta function. To determine the elastic proper- 
ties of the effective medium we adjust the medium spring 
constants a^ so that {bu) = 0, i.e. the lattice displace- 
ment in our spatially homogeneous effective medium ma- 
terial is identical to the average displacement in the spa- 
tially heterogeneous disordered material. 

Using this procedure we arrive at a spatially homoge- 
neous effective medium having spring constant am given 
by 



a 



if p > a* 



_ ) l-a* 

1 if p < a* 



(8) 



The contribution of network bending to the effective 
medium spring constant (which is proportional to the 
collective shear modulus) arises only through the effect 
of the bending modulus on a* in Eqs. (|3l4l5p . To de- 
termine how the shear modulus depends on the aver- 
age filament length we note the mean filament length 
is (i) = pi?(2 - p)/(l - p). We plot in Fig. [2] using the 
filled symbols the effective medium spring constant as a 
function of mean filament length measured in units of i?. 

We now consider the response of the network to a q- 
dependent strain as depicted in Fig. [TJ We modify both 
the bending modulus and spring constant of one filament 
spanning lattice sites so that: Km — > K'^am ct' 
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FIG. 2: (color online) The effective medium spring constant 
Qfm (filled symbols) and bending constant Km (open symbols) 
with the average filament length (L) (legend shows different 
values of k, with a = f). 



and compute the virtual force and torque needed to main- 
tain the position of site i in the middle of this triad of 
lattice sites (Fig.[T]). Using these forces and linearity, we 
compute the displacement of the i^^ site in response to 
the elastic constant subsitution made above. We find the 
displacements along the filament and perpendicu- 

lar to it (^£_l) are given by 

_ (ctm - o'){Uij + Uih).fij 

" 2{am/a* ~ a,n + a') 

where a* is defined in Eq. ^ and the analogous quantity 
b* is defined by 



i^Tr [D,{q)D-\q)] 



(10) 



using the same sum over wavevectors as in Eq. ([3]). Here, 
for semiflexible filaments on a triangular lattice interact- 
ing via cross-links that do not apply torques, the stretch- 
ing and bending modes are orthogonal p^ . 

We now average these displacements over the disorder, 
and, to find the effective medium elastic constants, we 
demand that the disorder-averaged displacements (^^||) 
and vanish. The probability distribution for a' is 

given by Eq. ([7]), but a nonzero value of the bending 
modulus at site i requires the presence of both filament 
segments on either side of that site so that 



P(k') = p-^Sin - k') + (1 - p'')S{n'). 



(11) 



Since we consider uncorrelated distributions of the 
bending and elastic constants we find the effective 
medium elastic constants am and Km by solving Eq. ^ 
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FIG. 3: (color online) The effective medium mechanical phase 
diagram spanned by (L) and The thick solid line marks 
the rigidity percolation transition where the material acquires 
a finite shear modulus. The dashed line shows the crossover 
from the non-affine to the affine regime. 



for ((5^11 ) = {Si±) = independently to arrive at 



a 



a p> a* 
ifp<a* 




(12) 
(13) 



Figure [2] shows the effective medium values of am 
(filled symbols) and Km (open symbols) as a function 
of (L) for different values of bending rigidity k, at fixed 
value of a = 1. The unit of length is the lattice constant 
R (set to unity) and the energy scale is arbitrary. 

There are three length scales in the system: (i) the 
average length of filaments (L), (ii) a length Z;, — [K/aY^'^ 
associated with the relative ease of filament stretching to 
bending, and (iii) the mean distance between cross links, 
which, to a good approximation is i?[14j|- We present 
a mechanical phase diagram of our system spanned by 
(L) and h in Figl3]that shows the regimes corresponding 
to zero and finite values of Km and am- Generically, 
for long enough filaments the system has finite collective 
extension and bending moduli. As the mean filament 
length is reduced the collective shear modulus vanishes 
at the rigidity percolation transition [13, [HI, ■ We also 
find a new rigid phase {am > 0) that has a vanishing 
bending modulus Km- We further note that the lower 
range in 4 corresponds to non-thermal systems in which 
the distance between crosslinks is large compared with 
the molecular scale (i.e., at low volume fraction), but 
that with thermal effects, one effectively goes to higher 

The collective elastic properties of the effective 
medium can be calculated from the stored elastic energy 
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FIG. 4: (color online) The effective medium shear Gcff/Go 
(filled symbols) and bending moduli K^s/Ko (open symbols) 
normalized by their respective values for a perfect network are 
plotted as a function of the mean filament length (L) divided 
by the nonaffinity length A = R/{R/lbY with z = 0.25. Data 
collapse is shown for four data sets differing in The legend 
shows different values of 



density £ under a given imposed network strain. For 
strain field of the form u = R'ycos{q ■ x)z x q the shear 
Geff and bending moduli ifeff can be extracted as the 
coefficients of the and q'^ terms of (f)/7^ where the 
angled brackets imply an angular average of the direction 
of q with respect to underlying lattice. Geff is propor- 
tional to Um alone while is a function of and 

In Fig.|4]we plot the effective medium shear and bend- 
ing moduli as a function of {L). Motivated by earlier 
work[5] on the A/NA transition we have rescaled (L) by 
A = R{R/lhY with z = 1/4. A comparison of Ges in 
this figure with am(oc Geff) from Fig. [5] demonstrates 
a remarkably accurate data collapse whose accuracy is 
enhanced as we move farther away from rigidity perco- 
lation. Moreover, we find that the same rescaling factor 



generates an equally accurate collapse of the KcS data. 

The collapse of our calculated elastic moduli with a 
single parameter, the lengthscale A, is in good accord 
with the numerical data collapse found in previous simu- 
lations Thus, the mean-field theory demonstrates all 
previously observed mechanical signatures of the A/NA 
cross-over seen there. The analytic results, however, 
suggest z « 1/4, while the prior numerics pointed to 
z = 1/3. The effective medium approach introduced here 
does not allow us to explore the spatial heterogeneities 
of the strain field under uniformly imposed shear, so it 
is impossible to explore the geometric interpretation of A 
with this technique. 

Although the efi^ective medium approach fails to ac- 
count for the correct spatial structure of the strain field 
in the disordered material, it does show an abrupt cross- 
over that appears mechanically identical to the A/NA 
cross-over. The cross-over is controlled by a single emer- 
gent length scale. A, which obeys a similar scaling re- 
lation to that found empirically from previous numeri- 
cal results. From these Geff plots (Fig. [3]) we have ex- 
tracted the A/NA cross-over from the location of the 
largest change in the slope of the curves. This A/NA 
boundary is plotted in Fig. [3] 

In conclusion, we used an effective medium theory to 
explore the mechanical properties of disordered filament 
networks. We find that this mean-field approach to the 
mechanics of such networks captures the mechanical as- 
pects of the A/NA cross-over including the identification 
of an emergent mesoscopic length scale A controlling the 
mechanics of the system. 
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